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Abstract 

This paper analyzes an abstract two-level algorithm for hybridizable discontinuous Galerkin 
(HDG) methods in a unified fashion. We use an extended version of the Xu-Zikatanov (X-Z) 
identity to derive a sharp estimate of the convergence rate of the algorithm, and show that 
the theoretical results also apply to weak Galerkin (WG) methods. The main features of our 
analysis are twofold: one is that we only need the minimal regularity of the model problem; the 
other is that we do not require the triangulations to be quasi-uniform. Numerical experiments 
are provided to confirm the theoretical results. 
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1 Introduction 

The Hybridizable Discontinuous Galerkin (HDG) framework, proposed in [17] (2009) for second 
order elliptic problems, provides a unifying strategy for hybridization of finite element methods. 
The unifying framework includes as particular cases hybridized versions of mixed methods 13 [3 
H5] . the continuous Galerkin (CG) method [03, and a wide class of hybridizable discontinuous 
Galerkin (HDG) methods. Here hybridization denotes the process to rewrite a finite element 
method as a hybrid version. It should be pointed that the Raviart-Thomas (RT) [32] and Brezzi- 
Douglas-Marini (BDM) mixed methods were first shown in [2; [7] to have equivalent hybridized 
versions, and an overview of some hybridization techniques was presented in M- In the so-called 
HDG methods following the HDG framework, the constraint of function continuity on the inter- 
clement boundaries is relaxed by introducing Lagrange multipliers defined on the the inter-element 
boundaries, thus allowing for piecewise-independent approximation to the potential or flux solution. 
By local elimination of the unknowns defined in the interior of elements, the HDG methods finally 
lead to symmetric and positive definite (SPD) systems where the unknowns are only the globally 
coupled degrees of freedom describing the Lagrange multipliers. We refer to mmm for the 
convergence analysis of several HDG methods for the second order elliptic problems. 

Closely related to the HDG framework is the weak Galerkin (WG) finite element method 
[34 1 (29] ]30j [31] pioneered by Wang and Ye 0J;. The WG method is designed by using a weakly 
defined gradient operator over functions with discontinuity, and then allows the use of totally 
discontinuous piecewise polynomials in the finite element procedure. By introducing the discrete 
weak gradient as an independent variable, as shown in E3, the WG method can be rewritten as 
some HDG version when the diffusion-dispersion tensor in the corresponding second order elliptic 
equation is a piecewise-constant matrix. 

It is well-known that the design of fast solvers is a key component to numerically solving 
partial differential equations. For the HDG methods as well as the WG methods, so far there are 
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only limited literature concerning this issue. In [22] (2009), Gopalakrishnan and Tang analyzed 
a V-cycle multigrid algorithm for two type of HDG methods for the Poisson problem with full 
elliptic regularity. By following the same idea, Cockburn et al. [T9] (2014) presented the first 
convergence study of a nonnested V-cycle multigrid algorithm for one type of HDG method for 
diffusion equations without full elliptic regularity. Chen et al. HU (2014) constructed two auxiliary 
space multigrid preconditioners for two types of WG methods for the diffusion equations. In [22] Li 
and Xie proposed a two-level algorithm for two types of WG methods without full elliptic regularity, 
and, in [25], they analyzed an optimal BPX preconditioner for a large class of nonstandard finite 
element methods for the diffusion equations, including the hybridized Raviart-Thomas and Brezzi- 
Douglas-Marini mixed element methods, the hybridized discontinuous Galerkin method, the Weak 
Galerkin method, and the nonconforming Crouzeix-Raviart element method. 

In this paper, we shall propose and analyze an abstract two-level algorithm for the SPD systems 
arising from the HDG methods for the following diffusion model: 

/ —div(aVu) = / in D, 

\ u = 0 on dfl, ' 

where fl C R d (d = 2, 3) is assumed to be a bounded polyhedral domain, the diffusion-dispersion 
tensor a £ [L°°(D)] dx<i is a SPD matrix and / £ L 2 (fl). In the two-level algorithm, the H 1 - 
conforming piecewise linear finite element space is used as the auxiliary space. The main tool of 
our analysis is an extended version of the Xu-Zikatanov (X-Z) identity [22] - The main features of 
our work are as follows: 

• We only need the minimal regularity of the model problem m in the sense that the regu¬ 
larity estimate 

Nl 1 +a,n < c 'n||/|| 0 _ ll n (1-2) 

holds with a £ [0,1], where Cn is a positive constant that only depends on D and a. Based 

on the convergence results of the two-level algorithm, Algorithm [I] (cf. Theorems 13. lllTTil) . 
it is easy to show that the multigrid methods which fall into the proposed two-level algorithm 
framework for the HDG methods all converge. We note that the analyses in m and m 
require full regularity (a = 1) and a £ (0.5,1], respectively. 

• We only assume the grids to be conforming and shape regular. Thus, the quasi-uniform 
condition, which is assumed in [221 HSl [HI O EH] > is not required in our analysis. Therefore, 
based on fast solvers for the auxiliary space, our analysis can be used to design fast solvers 
on adaptively refined grids and completely unstructured grids. 

• Our theoretic results also apply to the WG methods. 

The rest of this paper is organized as follows. Section 2] introduces notations, an extended ver¬ 
sion of X-Z identity and HDG methods. Section [3] describes and analyzes the two-level algorithm. 
Section [T] presents some applications of the algorithm to the HDG methods as well as to the WG 
methods. Section [5] reports some numerical results to verify the theoretic results. 


2 Preliminaries 

2.1 Notations 

For an arbitrary open set D C R d , we denote by H 1 (D) the Sobolev space of scalar functions on D 
whose derivatives up to order 1 are square integrable, with the norm H-^ D . The notation | ■ |i t £) 
denotes the semi-norm derived from the partial derivatives of order equal to 1. The space Hy(D) 
denotes the closure in H 1 (D) of the set of infinitely differentiable functions with compact supports 
in D. We use and (• ,-)od to denote the Ir-inner products on the square integrable func¬ 

tion spaces L 2 (D) and L 2 (dD), respectively, with ||-|| D and ||-|| 9£) representing the corresponding 
induced L 2 -norms. Let Pk(D) denotes the set of polynomials of degree < k defined on D. 
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Let 7ft, be a conforming and shape regular triangulation of f2. For each T £ Th, hr denotes 


the diameter of T with h := ma xtgTh The regularity parameter of 7ft is defined by p := 

max-reTh ^t/|T|, where \T\ is the d-dimensional Lebesgue measure of T. Let J -ft denote the set of 
all faces of 7ft. 

We define the mesh-dependent inner product (•, -)h and the corresponding norm [|-||, as follows: 
for any X,p£ L 2 (JF h ), 

(X,y)h ■= ^2 / 'Vb 

TeT h JdT 

Uplift := (A,A)J. 

(2.1) 

We also need the following notations: 



Uplift,ar := Hollar 

, VA e L 2 (dT), 

(2.2) 

\k\h '■= ( ^2 H ft.ar) 2 

, Vp£L 2 (T h ), 

(2.3) 


T£T h 


where 

H l,dT '■= h T X IlM - m T {^)\\ 2 aT , rn T (n) ~ f dT p, 

and \dT\ denotes the (d-l)-dimensional Hausdorff measure of dT. 

Throughout this paper, x < y (x > y) means x ^ Cy (x ^ Cy), where C denotes a positive 
constant that only depends on d , k , f2, the regularity parameter p, and the coefficient matrix a. 
The notation x ~ y abbreviates x <y <x. 


2.2 Extended version of X-Z identity 


We start by introducing some abstract notations. Let V be a finite dimensional Hilbert space 
equipped with inner product (•,•) and its induced norm ||-||. Suppose A : V —> V is a linear 
operator which is SPD with respect to (■, •), then (-, •) A '■= {A-, •) also defines an inner product on 
V and we use ||-|| A to denote the corresponding norm. Let B : V —> V be a linear operator with 
norm 


\\B\\ a := sup 

vGV 


\\Bv\\ a 

ll«L ' 


Suppose Vo, Vl,V 2 ,- • • ,V/v are finite dimensional Hilbert spaces equipped with inner products 
(•, -)o, (•, -)i, • • •, (•, -)jv respectively. Let Ii : V) —> Vo (i = 1, 2, • • • , TV) be linear injective operators 
such that 


V 0 = hVx + l 2 V 2 + --- I N V N . 


Naturally, the adjoint operator 7| of Ii is defined by 


(IiVo,Vi)i = (v 0 ,IiVi) 0 for all Vi £ V and v 0 £ V 0 . 


Let Ao : Vo —t Vo be SPD with respect to (•,-)o and define A, : Vi —► Vi (i = 1,2, ,7V) 

by A, = I*Aoh- Since Ii is injective, A t is SPD with respect to (•,•)*• For each i, suppose 
IZi : Vi Vi is a good approximation of Ay 1 and define the symmetrization of 77,, by 

IZi = IZ* + IZi ~ IZ.]A, IZ ,. (2-4) 


Then we define the operator Bq : Vq —> Vq as follows: 


For any given b 

£ Vq , B 0 b := v 2N with v 2N defined below: 

v° := 0. 


for 7 = 1,2,- 

• ,N 

v l := v 1 - 1 

+ Ii1ZiI*(b - Aqv*- 1 )] 

end 


for i = TV, TV 

— 11 

v 2N-i+l 

= V 2N ~* + IiIZ\Il{b - A 0 v 2N ~ i y, 

end 



Finally, following [3H1 HI El HU, we are ready to present the following extended version of X-Z 
identity. 
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Theorem 2.1. Suppose IZi is such that ||/ — 7^.4; ||^ < 1 for i = 1, 2, • • • ,N. Then it holds 


\\I-BoMa 0 



(2.5) 


where 



3>i 


Hi 


( 2 . 6 ) 


Proof. The desired result can be obtained by following a similar routine to the proof of (Theorem 
4.1, [53]), which is a trivial modification of the new proof [5] of the X-Z identity. □ 


2.3 HDG framework 

We give a brief description of the HDG framework; One may refer to m for more details. For 
any T £ Th, let V(T ) C L 2 (T) and W{T) C [ L 2 (T)] d be finite dimensional spaces, be a 
nonnegative penalty function defined on dT , and P|( : L 2 (dT) —> M(dT) be the standard L 2 - 
orthogonal projection operator with 

M(dT) = {pL £ L 2 {dT) : p\p £ Pfc(F), for each face F of T}. 


Introduce the finite dimensional spaces 

V h := {v £ L 2 (Cl) : v h \ T £ V(T), VT £ 7^}, (2.7) 

W h := {t £ [L 2 (n)] d : t^| t £ W(T), VT £ T4, (2.8) 

Mh := {At* £ L 2 (F h ) : /z^f £ Pfc(F), VP £ and /x/j|ao = 0}. (2.9) 

The general framework of HDG methods for the problem (11.11) reads as follows ([12): Seek 
(uh, Xh,er h ) €V h x M h x W h such that 

(' Ctr h , Th) + (u h , divhTh) - T. (A h , T h • n) dT = 0 X/t h £ W h , (2.10a) 

TdT h 

-(v h ,div h cr h ) + ^2 (a T (PTUh ~ X h ),v h )oT = (f,v h ) Vv h £ 14, (2.10b) 

Ter fe 

^ (c/i • n - ariPrUh — Xh), Th)dr = 0 V/z^ £ M h . (2.10c) 

TGTh 


where c = a 1 , and divh is the broken div operator defined by divhTh |t := div(Th\T) for any 
T h G W h ,T £ %■ 

Introduce the following local problem: for any A £ L 2 (dT), seek (u x , cr x ) £ V(T) x W(T) such 
that 


(c<t a ,t) t + {u\,divr) T = (A ,T-n) dT Vt £ TT(T), (2.11a) 

— (i>, divcr\) T + (arPrUx, v)ot = {arX,v)dT Vi> £ V(T). (2.11b) 

Let ah {•, •) : M/, x M/j -^Ibea bilinear form associated with the above local problem, defined by 

a h (X h , fi h ) := ^2 (c(T Xh ,(T flh )T + J2 MP$u Xh - X h ),P$u, h - p h ) dT . (2.12) 
TgTh TeTh 

Then the HDG model (12.101) is equivalent to the following reduced system [IT]: seek A h G Mh such 
that 

Q J h{Xhi Th') = (f l’U’Uh) o, V/Zfo G Mh - (2.13) 

We note that once the Lagrangian multiplier approximation A h is resolved, the numerical flux crh 
and the potential approximation Uh can be obtained in an element-by-element fashion by (12.111) . 
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3 Two-level algorithm 

We recall that the triangulation 7 ~h is assumed to be conforming and shape regular. In addition, 
we assume the regularity estimate CL2D holds with a £ [0,1]. 

For the sake of convenience, in the rest of this paper we shall use the notation (•, •) to abbreviate 
the L 2 -inner product (•, -)q. 

3.1 Algorithm description 

At first, we introduce the it^^-conforming piecewise linear finite element space 

V h := {v h £ H^Q) : v h \ T £ Pi(T)}. (3.1) 

We then dehne the prolongation operator A : 14 —» Mh and its adjoint operator I f h : Mh —► 14 
respectively by 

(ThVh , t^h) h — {vh,t^h)h 'tfvh £ 14 ; ( 3 . 2 ) 

{J-hVhi ^h) ~ {t^h, ThVh)h £ Mh, ( 3 . 3 ) 

and define the operators Ah ■ Mh —> Mh and Ah ■ 14 —> 14 respectively by 

(AhXh,Hh)h = ah{^h,^h) (3.4) 

(■ A h u h ,v h ) = (aVtifc, Vuh) \/u h ,v h £V h . (3.5) 

Let IZh ■ Mh —> Mh and 7 Zh : 14 —> 14 be good approximations of A jj" 1 and Ah respectively, 

— t 

with TZ h and IZh satisfying 

Vh)h ~ {Xh,TZh^h)h *TXh, t^h £ Mh, 

(7 Z h u h ,v h ) = ( u h ,TZ h v h ) \/u h ,v h £ V h - 

Finally we define the operator Bh : Mh —> Mh as follows: 

For any rjh £ Mh, BhJ]h ■= Hh with ^ dehned below: 

1. [ := TZ h r]h', 

2- Vh := Vh A IiiTthlhirih — Ah^\); 

3 - Vh : = vl + ihR-h Ihivh - A h nl); 

4 - vi ■= V 3 h + K(Vh - AhUp. _ 

In view of the operators Ah and Bh, we present the following two-level algorithm for the system 


E 

3: 


Algorithm 1. Let bh £ Mh be given. We solve the equation AhXh = bh below: 

A£ = 0, 

for j = 1, 2, • • • 

Xi^Xi^+Bhibh-AhXr 1 ); 

end 


3.2 Main results 

We first introduce the following symmetrizations of TZh and 7 Zh'- 

IZh '■= TZ], + TZh — IZhAhlZh, 

TZh := TZh + TZh — TZh AhTZh, 

TZh ■= TZh + TZh —TZh AhTZh, 

where _ 

Ah := IhAhlh- 

Then we present some assumptions below. 


(3.6) 

(3.7) 

(3.8) 

(3.9) 
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Assumption I. For any A h £ Mh, it holds 

j9 


1 

2 

C2(T \ h 

+ 

T 


Assumption II. It holds 


where 


M h := 


1 -M> 


I - IZhAh 


a^(P T U\ h A h) ~ \^ h \h,dT > VT € Th- 

(1 + Mh)Mh +Afh 


1 + M h 


■Ah 


Mh ■= 


< 1, 


I - Ah Ah 


A h 


Assumption III. Let IZh '■ Mh —> Mh be SPD with respect to (•, •)h such that 

0 < cr(TZ h Ah) < w, 

_1 

l|A h ||^- 1 < (EreT.^II^TY, VA h £ M h , 


(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 


where a(IZhAh) denotes the set of all eigenvalues of IZhAh, and uj is a constant with 0 < ui < 2. 

Remark 3.1. Obviously, Assumption U implies 

ah(Ah, M) ~ \Ah\h , VA h £ M h , 


and hence Ah is SPD with respect to (-,-)h- Since Ih ■ 14 —>• Mh is injective, it is easy to verify 
that Ah is SPD with respect to (-, •). What’s more, by a simple estimate | Aft, dT < hf, 1 ||Aft|| 

h,OT’ 

Assumption [I] implies that the largest eigenvalue of Ah satisfies A max (Ah) ^ h~ 2 under the 
condition that 7 h is quasi-uniform. 

Remark 3.2. In Assumption HH when Mh is given, the condition (13.111) requires that A ih is 

sufficiently small, i.e. the operator IZh is a good-enough approximation of Ah ■ Fortunately, for 
the H 1 -conforming linear element approximation Ah, the research of the choice of IZh is mature. 

As will be shown in Section^ for some applications, it holds Mh = 0 or Mh fz, h. In the former 
case, (13.111) is reduced to the constraint 

M h < 1. (3.15) 

In the latter case, h should be also small enough to ensure (13.111) . We note that Assumption HT1 
requires implicitly the constraint (13.1511 . 

Remark 3.3. It is evident that the condition (13.131) in Assumption Hill implies that ||/ — IZhAh ||^ fc < 
1, which means 

||I - AZhAh\\ Ah = ||/ — 7ZhA h \\ 2 Ah < 1. 

Suppose Assumption|T]*s true. If we choose the Richardson iteration as IZh, i-e. IZh = j — 1 (A h ) ^’ 
then (13.131) holds with u> = 1, while (13.141) holds only in the case that Th is quasi-uniform. However, 
if we set IZh to be the symmetric Gauss-Seidel iteration, then (13.131) holds with ui = 1, and (13.141) 
holds as long as Th is conforming and shape regular. We refer to Appendix A for a concise 
analysis of the symmetric Gauss-Seidel iteration. 

We state the main results in two theorems below. 


Theorem 3.1. Under Assumptions m um it holds 

\\{I-B h A h )\\ Ak = \-^, 


where 


K<l+M h + 


(1 + Mh )/(2 — LO) 


1 1+Wh 

1 ~ TWW 


1 2 ' 


(1 + Mh)Mh + Mh 


(3.16) 

(3.17) 
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Theorem 3.2. Let TZh be one sweep of Gauss-Seidel iteration. Then, Under Assumptions mm 
the relation (13.1611 holds with u> = 1. 

We shall prove these two theorems in Section 13.31 

Remark 3.4. Since we only assume Th to be conforming and shape regular, it’s important that 
Theorems \3. 1 f3f2\ hold on non-quasi-uniform grids, as long as there is a proper choice of IZh for the 
H 1 -conforming linear element approximation. We refer the reader to MWiUWMUW 

for the construction of IZh on adaptive grids, and to EUBEauSP for the construction on completely 
unstructured grids. 

3.3 Convergence analysis 

3.3.1 Proof of Theorem 13.11 

Lemma 3.1. Under Assumptions HUTTl it holds 


I — IZh Ah 


m < 1 > 

Ah 


--i ^ 


1 + Afh 


Tlh 


1 - 


l±Wh 

TASK 


1 2 


(1 + J\fh)A4h + Afh 


IMfo, Vv h &Vh. 


Proof. Since Ah Ah is symmetric with respect to (■, , we have 

(1 ~Af h ) IKII^ < |K|||= < (1 +A fh) |K||~ , Vv h e V h . 

For any linear operator Sh : 14 —> Vh, it holds 


(3.18) 

(3.19) 


(3.20) 


\\S h Vh\\~ 

ll^ll- = sup ■ 

Ah v h ev h 




1 + A fh 


1 -Arhv h Jv h IK 11^ 

S h \' 


IIKKI 


A h 


1 + Afh 


l-Af h 


"Ah 


(3.21) 


Then, from 


it follows 


/ - TZhAh 

< 

I — TZhAh 

_ + 


A h 


Ah 


TZhAh {I - Ah Ah) 


^ Aih A Afh 


TZhAh 


(1 + Afh)A4h + Afh 


A h 


I - TZhAh 


I — TZhAh 


Ah 



€ 


€ 


1 + Afh 

1-A f h 

1 + Afh 

1-A f h 


A h 


I - TZhAh 


Ah 

(1 + Afh)A4h + Afh 


(by (13.211) 1 
2 


(3.22) 
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which immediately implies (13.1811 . 

On the other hand, by the definition of IZh , we can get 


I - TZhAh 


— 1 A min (7 ZhAh), 


Ah 


where \ m in{TZhAh ) denotes the smallest eigenvalue of TZhAh- The above relation, together with 
the fact that, due to (13.181) . IZh is SPD with respect to (-, •), yields 


A min(TZhAh ) 


1 


1 - 

/ — TZhAh 




A h 


IKII 


(3.23) 


Ah 


Finally, the desired inequality (13.191) follows immediately from (13.201) and (13.221) . This completes 
the proof. □ 


Lemma 3.2. Under Assumptions Ul lIIIL the relation (13.101) holds with 


K = sup 


IP h\\ Ah =i^ h+IhVh = Xh 


inf ^ || fi h + TZhAhlhVh 11^-1 + IK II—- 1 ■ 


Hh 


(3.24) 


Proof. The conclusion follows from the space decomposition Mh = Mh + IhYh and the extended 


version of X-Z identity (12.51) . 

Lemma 3.3. Under Assumptions Ul lIIIl it holds 

\\TZ h AhIhVh\\^-i < max j 1, 0 U _ [>(1 + A f h ) IKIIjb , Vv h G 14. 


2 — uj 


□ 


(3.25) 


Proof. Denote Sh := TZhAh- Noting that 


A4 = 272,/j — TZhAhTZh = (2<54 — Sf,)A h 1 . 


we have 


\\TZhAhIhVh\\f^-i = (714 TZhAhIhVh,TZhAhIhVh)h 


= ( Ah{2Sh - Sh) 1 ShIhVh,TZhAhIhVh)h 
= (■ Sh{2Sh - Sl)~ 1 ShIhVh,IhVh)A h - 


(3.26) 


Since <54 is SPD with respect to (•, -)^ h and the inequality 


t(2t — t z ) 1 t ^ max < 1 


w 


2 — oj 


t G (0, Cel] 


holds, the relation (13.261) . together with (13.201) . immediately yields the desired estimate (13.251) . □ 

From Lemmas I3.H3.3I and Assumption IIIII we obtain immediately the lemma below. 


Lemma 3.4. Under Assumptions Ul lIIIl the relation (13.161) holds with 

(1 + Mh )/(2 — to) 


K< sup inf h T 2 \\nh\\l dT + 


\\\ h \\ Ah =l^ +Ihvh=xh T&Th 


1 1 +A 5 ,, 

1 ^ 


1 2 


(1 +Mh)Mh +Mh 


■Ah 


(3.27) 
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To further derive (13.171) . we introduce the operator Ph : Mh —> Vh with 
P h X h (x) = < 


E rn T (\h) 

, if ai is an interior node of Th, 


E i 

Tei^a, 

0, if x G dfl 

for any A h G Mh, where denotes the set {T G Th ■ x is a vertex of T}. We have the following 
important estimates for Ph- 

Lemma 3.5. For any A h G Mh, it holds 

\P h Xh\i,n 


\\(I ~ hP h )X h \\l 


dT 


< \Xh\ h , 

(3.28) 

< |A,| h . 

(3.29) 


Proof. For each T £ Th, we denote wr := {T' G Th : T’ and T share a vertex} and use M(T) to 
denote the set of all vertexes of T. Assume all vertexes of T are interior nodes of Th, then it holds 


\\PhXh — m.T(Xh)\\g T 




< h T ExejV(T') \PhXh{x) m T {X h )\ 2 

scejV(T)E T u T a €u>„ \m Tl {\h) ~ m T2 {\h)Y 

Ti,T 2 share a same face 


~ Exejvcr) E Ti^eua, WmrAXh) - mT 2 {Xh)\\ dTl ndT 2 

Ti,T 2 share a same face 

~ Et'Gwt h T ' I Xh\h,dT> ' 


(3.30) 


Similarly, we can show by a trivial modification that (13.301) also holds in the case that there is a 
vertex of T that belongs to dQ. As a result, the estimate (13.281) follows from 

\PhXh\i,T = \PhXh — rriT^Xh) If x 

< hf 2 \\PhXh — mT(Xh)\\x (by inverse estimate) 

< hf 1 \\PhXh — mT{Xh)\\ 2 gT 

< Et'Gwt I ^h\h,dT’ ■ (by ( | 3 . 30 p ) 

On the other hand, from 

hr \\IhPhXh — rriT{Xh)\\ 2 dT < hr \\PhXh — wt(A/ 1 )||| t 

Er'e^r h T- \ X h\l, aT ' (by (|3.30p) 


it follows 


hx ||(/ — IhPh)Xh\\ 2 dT < hr ||Ah — TOr(Aft)||g T + hr \\IhPhXh — mr(XT " 2 


I9T ~ |\''h '"‘1 V'nJMdT 

< ^ hx> \Xh\h,dr’ > 
T'Gcjt 

which indicates (13.291) immediately. 


dT 


□ 


Remark 3.5. Although similar estimates were presented in \2‘J[ 1 25/ for quasi-uniform grids, the 
estimates in Lemma \3.5\ are sharper in the sense that Th here is not assumed to be quasi-uniform. 

Finally, we are in a position to prove Theorem 13.II 
Proof of Theorem l3.ll For any A h G Mh, set ph ■= Xh~ IhPhXh and Vh := PhXh ■ Using Lemma 
13.51 we have 

(1 + -A/*.)/(2 — u) 2 

-T2 IKII A h 


Z ,l 

T&T h 


-2 

T 


2 

h,8T 


1 - 


l+Afh 


(1 + A fh)A4h + Mh 


< < 1 


(1 + Mh )/(2 — ui) 


1 - 


l+ACh 

U7v7 


(1 + Mh)Mh + A/} 


IIA/J 


> 


9 

























which implies 


s U p inf Y^ h T 11^1 \ h ,ar + — 

\\^U h =^ h+IhVh = Xh ^- h 

1 Wa/I 


(1 + A// l )/(2 — w) 


(1 + Hh)M-h + .A/h 


2_ 


<1 + 


(1 + Mh )/(2 — w) 


1 - 


l+A/fe 

T^vTT 


(1 + Mh)Mh + A/h, 


2 • 


Then Theorem 13.11 follows from Lemma 13.41 immediately. 


3.3.2 Proof of Theorem 13.21 

Let {ip : i = 1,2,..., A/ - } be the standard nodal basis for Mh . We have the following space 
decomposition: 

Mh = span{r/i} + span{? 72 } + ... + span-f^} + I h V h . 

Define Pi : M h ->• span{? 7 j} by (PiX h , rji) Ah = (A h , rji) Ah for i = 1,2,..., N. Then, by the extended 
version of X-Z identity (12.51) . we have the following lemma. 

Lemma 3.6. Under the conditions of Theorem 1 3. 2\ the relation (13.161) holds with 


N 


K = sup inf y 

\\\ h \\ Ah =ii:!Li iH+ihv h =\ h ^ 


Pi + Pi l ^2 Pj + IhVh ) 
A j>i / 


Tth 


A h 


Lemma 3.7. Under Assumptions HUTT1 it holds 


N 

E 


Pi 


Pi ( ^2 Pj + IhVh 


j>i 


< 


E h T 2 \\Ph\\ h ^ dT + (1 + A fh) ||w/i|| 


2_ 

Ah 


Ah 


TeTh 


(3.31) 


(3.32) 


for any v h £ V h and m £ span{rji} (i = 1,2,..., N) with p h = Pi- 

Proof. Define E, := [T : there exists one face F of T such that ry \p ^ 0} and w, := Ures*? 1 . 
Apparently, for any given 1 ^ io ^ AT, there are at most J of {w^}, (w,. : j = 1,2,..., J}, such 
that uji 0 D u>i j ^ (j) (j = 1,2,..., J), where J only depends on the dimension number d and the 
shape regularity parameter p. 

It is easy to verify 


N 


y\\phvh\\ 2 Ah <\\i h v h \\ 2 Ah 


i=l 


(3.33) 


<(1+Mh)\\v h \\~ h (by») 


and 


N 


E 


Pi + Pi E Tj 

2 

II 

Pi+Pi y Pj 

j>i 

Ah 

i =1 



N 


< 


E ii^ 




Ah 


N 


< 


E I A 4 * I h- (by Assumption [T]) 


i— 1 


Then, from 


\p\h,dT ~ M\h,dT > ^ M(dT) 


(3.34) 
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it follows 


N 

E 

i =1 


fj'i t p% y \ f-^j 

j>i 


where we have used the estimate 


N 


N 


< 


< 


E E ^ t 2 IlMtll/j.ar 

*= i TeTh 


JV 


E E 11^*11 h,dT 

TeTh i=l 


< ^ /i T 2 ||/X 

TeTh 


h \\h,dT > 


E ~ IlMfclU.dT i 


(3.35) 


which can be proved through standard scaling arguments. Consequently, the desired estimate 
(13.321) follows immediately from (13.331) and (13.351) . □ 

By Lemmas 13.6113.71 and (13.191) , we immediately obtain the lemma below. 

Lemma 3.8. Under the conditions of Theorem \3. c A the relation (13.161) holds with 

1 + Nh 


K < sup inf h T 2 

\M\A h =^ h+IhV h = X h r^- h 


2 

h,dT ' 


1 - 


1 +A /"h 


(1 + + Mb, 


•Ah 


. (3.36) 


Finally, the rest of the proof of Theorem 13.21 goes exactly the same way as that of Theorem 13.1 


4 Applications 

This section is devoted to some applications of the algorithm analysis in Section l5~21 to some existing 
HDG methods as well as WG methods. 

In the two-level algorithm, Algorithm [T) described in Section 13.11 we set the operator Rh 
to be the symmetric Gauss-Seidel iteration or one sweep of Gauss-Seidel iteration. As shown in 
Remark 3.3 and Appendix A, the symmetric Gauss-Seidel iteration always satisfies Assump¬ 
tion [TTT1 Thus, according to Theorems 1.3. llTt~2l we only need to verify Assumptions HUTT1 for the 
corresponding methods. 

We consider the following four types of HDG methods: For any T £ Th, k > 0, 

Type 1. V(T) = Pk(T), W(T) = [Pk(T)] d + Pk(T)x and or = 0. The corresponding HDG 
scheme (12.101) is the hybridized RT mixed element method ([2]). 

Type 2. V(T) = Pk~i(T), W(T ) = [Pk(T)] d (k > 1) and ax = 0. The corresponding HDG 
method is the hybridized BDM mixed element method ([7]). 

Type 3. V(T) = Pk(T), W(T) = [Pk(T)] d and ax = 0(1). The corresponding HDG method 
was proposed in m and analyzed in [18.. For the sake of simplicity, we assume for this HDG 
method that ax is constant on dT but it may take different values for different elements T. 

Type 4. V(T) = Pk+\(T), W(T) = [Pk{T)] d and ax = 0(hf}). The corresponding HDG 
method was analyzed in |24j 

For these HDG methods. Assumption [T| has been verified in mum for Types 1-3 methods 
and in [23] for Type 4 method. Then it suffices to verify Assumption [H] 

For the diffusion-dispersion tensor a, we consider two cases: piecewise constant coefficients and 
variable coefficients. 
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4.1 Piecewise constant coefficients 

In this subsection, we assume a to be a piecewise constant matrix, and, without lose of generality, 
we just take a to be the identity matrix, since the analysis is the same as that of the former case. 

Let w £ -Pi(T) and set A = P^w in the local problem (12.111) . For Types 1-2 HDG methods, 
it is trivial that 

cr\ = "Vw. 

For Type 3 (k ^ 1) and Type 4 HDG methods, we can easily obtain 

= A, er A = Vw. 

Thus, by the definitions (I3.2I) - (I3.5I) and (13.91) . for all the mentioned cases above we easily have 


■Ah = Ah , 


(4.1) 


I - Ah Ah 


(cf. (13.121) 1 and Remark 13.21 indi- 


A h 


which, together with the definition TV), := 
cates the following conclusion. 

Proposition 4.1. For Types 1-2, Type 3 (7c ^ 1) and Type 4 HDG methods, it holds 

A f h = 0, (4.2) 

which implies that any choice of IZh satisfying (13.151) ensures Assumption [TT] to hold. 

For Type 3 HDG method in the case k = 0, we have the following result. 

Proposition 4.2. For Type 3 (k = 0) HDG method, it holds 

A r h < h, (4.3) 

which implies that sufficiently small mesh size h can ensure Assumption [IT] to hold if IZh satisfies 


(13.151) . i.e. Mh < 1, with Mh '■= 


I — TZhAh 


Ah 


being independent of h. 


Proof. For any w £ Pi(T), set A = P^w in the local problem (12.111) . then it holds 


er a = Vro, u\ = mriw) = 


1 


1^1 JdT 

Consider an auxiliary problem as follows: for any Uh £ 14., seek Vh £ V), such that 


(4.4) 


(■ U h ,w h )~ = (v h ,Wh)jii £ Vh. 

A h 


By (14.41) . we easily obtain 


(Vu/,, Vw h ) + ^ (aT(rn T (u h ) - P^u h ),m T {w h ) - P^w h ) dT = (V®/,, Vw h ), Vw h £ 14, 
TeT h 

and it follows 

(V(u h - v h ), Vw h ) = - ^ (a T (m T {u h ) - P^u h ),m T (w h ) - P^w h ) d T, Vw ft £ 14. 

Since 


Ten 


||« - m T (v)\\ dT < hf || Vu|| T , Mv £ H\T), 


(4.6) 
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(by (gH)) 


taking Wh = Uh — Vh in (14.51) we have 




< 

< 

< 


D a T || mr{uh) - u h \\ aT || m T (w h ) - w h \\ aT 

TeTh 

ot T h T \\Vu h \\ T \\Vw h \\ T 

T&T h 

max a T h T \\u h \\^- h \\u h - , 


which leads to 


i.e. 


||Uh - Vh IIjar < max a T h T ||u/,|| 


TdT h 


"A h 


(4.7) 


(/ — Ah Ah)uh 


< max a T h T ||u/,||^ , Vu h e V h 


A h 


TdT h 


which, by recalling a-r = 0(1), yields (14.31) immediately. 


□ 


Remark 4.1. By Theorems \3. llTOl it is easy to derive the convergence rate (independent of mesh 
size) of a V-cycle HDG multigrid in 122f . where full elliptic regularity (Q was assumed to be convex) 
was required. However, our analysis does not require full regularity. 


Remark 4.2. From Theorem \3. 1\ and Proposition \f . 2j in order to the convergence of Algorithm 
[lj we have to require h to be small enough. This is in agreement with the theoretical result in fT9j. 

Suppose IZh satisfies (13.151) . We summarize this subsection as follows: 


• For Type 1-2, Type 3 (k ^ 1) and Type 4 HDG methods, both Theorem 13. 1 1 and Theorem 
13.21 hold with 

K < 1 + — 


1 


1 - 


I — 7 ZhAh 


A h 


• For Type 3 (k = 0) HDG method, the mesh size h should be sufficiently small to ensure the 
convergence of Algorithm [T] 


4.2 Variable coefficients 

In this subsection, we assume a £ [W 1,oc (Th)] dxd , where W 1,00 (Th) '■= {a € : Vu|t £ 

[L°°(T)] d , VT £ Th}- In the analysis below, we only consider Types 1-2 and Type 3 (k ^ 1) 
HDG methods, since by the technique used here, it is easy to derive similar results for other HDG 
methods. Following the same routines as in Section I4~T1 Icf. Propositions 14.1114.211 . we only need to 


estimate the number Afh = 


I - Ah Ah 


Ah 


Lemma 4.1. For Types 1-2 HDG methods, it holds 


Af h < h. (4.8) 

Proof. For any w £ P\{T), set A = P^w in the local problem (12.111) . Then it is easy to show 

M T <|A| MT ~||v™|| T . (4.9) 

On the other hand, by (12.111) we also have 

(ccr\ — Vw, ct\ - cT 1 Vw)t = 0 

with c := -ryj ■ f T c. Thus it holds 

(ccr A — Vw, cr\ — c~ 1 'Vw)t = ((c — c)cr\,(T\ — cT 1 Vw)t + (ccr A — Vro, cr\ — c~ 1 'Vw)t 

= ((c - c)(T\, CT\ - cT 1 Vw)t 
< h T ||er A || T ||er A - c _1 Vw;|| T , 
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where in the last we have used the standard estimate 

l|c ~ c|| £00(7^ < hr II Vc|| L <x,( T j 

\\cct x - Vw|| T < h T ||<x A || T , 


Hence it follows 
which implies 


||cer A - Vw|| T < ||(c — c)er A || T + ||cer A - Vic|| T 

^ hr Ho-aIIt- 

This estimate, together with (HU), yields 

||ccr A - Vw|| T < h T || Vw|| T . 


(4.10) 


(4.11) 


(4.12) 


Finally, for any Vh £ 14, taking w = vh\t in (14.91) and (14.121) with A = IhVh\r, from the definitions 
(I3.2M3.5I) and (13.91) . it follows 


{{A h - Ah)v h ,v h ) 


Ik/^JIc- ll aV ^llc 


4 (Ik/. 


h v h IIc 


+ IkV^IL) \\(T IhVh - aV%|| f 


;$ h bh\\ 2 A h > 

which gives the desired estimate (14.81) . 

Remark 4.3. For Type 1 HDG method (k = 0), if we redefine Ah as 

{A h u h ,v h ) = {c~ 1 'Vu h ,'Vv h ), Vu h ,v h £ 14, 

then it holds Ah = Ah and Afh = 0. This is a trivial modification of W- 

Next we consider Type 3 HDG method. 

Theorem 4.1. For Type 3 HDG method (k 4 1/. the estimate (14.81) holds. 

Proof. Let w £ P\{T) and set A = P®w in the local problem (12.111) . It is easy to obtain 

(cer A — Vie, t)t + {u\ — w, divT)r = 0, Vt £ W (T), 
-(v,diva\) T + (a T {u\ - X),v)dT = 0, \/v £ V{T). 

Taking t = t r A — c _1 Vw, v = u\ — w and adding (I4.13al) and (I4.13bl) . we have 

(ccr A - Vro, <x A - c~ 1 Vw) T + (a T {u\ - A),u A - \} dT = 0. 

This relation yields 

(cer A - Vu>, <x A - c~ 1 'Vw) T + (a T {u\ - A), u\ - A ) dT 
={{c - c)er A , <t\ - c~ 1 Vw) T 
<h T ||<t a || t ||< t a - c _1 Vic|| T , 

which implies 

1 

||cer A - Vw\\ T + afi ||u A - A|| 9T < h T ||<t a || t ■ 


Hence it follows 


□ 


(4.13a) 

(4.13b) 


(4.14) 


|cer A - Vin|| T + afi ||w A - A|| ar < h T ||er A || T , 


(4.15) 
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which, together with flUt , shows 


II ca x - Vio|| T + af, ||ma - A|| ar < h T || Vw|| T . (4.16) 

Finally, for any Vh £ 14, taking w = Vh\r in ( 021 ) and (14.1611 with A = IhVh |t, from the definitions 
m-m and it follows 


{{A h - A h )v h ,v h ) 


<h 


ik.. lie + £ a T || Ul h v h - IhVl i|| 
T&T h 

A h ’ 


2 

dT 


\\aVv h \\l 


which implies (14.81) . □ 

Remark 4.4. We note that our analysis only requires the regularity estimate CE2D with a £ [0,1], 
while the analysis in m requires a £ (0.5,1]. 

Similar to Section [4.11 we summarize this subsection as follows: 

• When the tensor a is not piecewise constant but piecewise smooth, the convergence of the 
two-level algorithm, Algorithm [T] for the HDG methods can still be obtained, as long as 
the mesh size h is small enough. 

4.3 Application to weak Galerkin methods 

In this subsection, we shall show our analysis can also be extended to the WG methods. Unless 
otherwise specified, we adopt the notations introduced in section [2] 

Following m , we introduce the weak gradient operators as follows. For any T £ T h . define 
V*, : L 2 {T) -+ W{T) by 

(V>, q)T = -(«, divq ) T , Vq £ W{T),Vv £ L 2 {T), (4.17) 

and V b w : L 2 {dT) -> W{T) by 

(Vt A, q) T = (A, q ■ n) dT , Vq € W(T), VA £ L 2 {8T), (4.18) 

where n denotes the unit outward normal vector to dT. 

The WG framework for the model problem (hud reads as follows ([Mj): seek {uh, A h) £ 14 x Mh 
such that 

y. {a{V l w u h + V b w X h ), V l w v h T V b w iu, h ) T + S h {{u h , \ h ), {v h , Uh)) = {f,v h ), V(v h ,fi h ) £ V h x M h , 

TGT h 

(4.19) 

where 

Sh{{uh, A h, {Vh, Uh)) '■= {otT{PT u h — A h),Pr v h — l^h)dT• (4.20) 

T&Th 

Denote <T/ t := + V^A/j, then the WG model (14.191) is equivalent to the following HDG-like 

scheme: seek {uh, Xh, <Th) £ Vh x Mh x Wh such that 

{erh, r h ) + {u h , divT h ) - ^ (A h , r h ■ n) dT 
T&T h 

-{v h ,div h {P^{acrh))) + ^2 (<*t(-PtKi “ A h),v h ))dT 

T&T h 

yy (P™{acT h ) ■ n - a T {PrUh ~ A h),Vh)dT 
TeT h 

where P™ : [L 2 {£l)] d —> Wh denotes the standard L 2 -orthogonal projection operator. 


= 0, Vr h £ W h , (4.21a) 
= {f,Vh),Vv h £\ 4, (4.21b) 

= 0, V/i h £ M h , (4.21c) 
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We define the local problem as follows: for any A £ L 2 (dT), seek (m“ 9 ,<t™ 9 ) £ V(T) x W(T) 
such that 



( <t a 9 > t )t + (u™ 9 ,divr) T = (X,t ■ n) dT , Vr £ W{T), 
-{v,div(P™(acr™ 9 ))T + {a T PTU™ 9 ,v) aT = {a T X, v) dT , Vu £ V{T), 

(4.22a) 

(4.22b) 

where P™ : [L 2 (T)] d —> W(T ) denotes the local L 2 -orthogonal projection operator. 
Similar to Theorem 2.1 in [T7], the following proposition holds. 


ProDosition 4.3. Suvvose (uh.Xh) £ 14, x Mh solves the WG model (14.191). then 
obtained by solving the system 

A h can be 


a™ 9 {X h ,Hh) = (J>™ 9 ), VfJ, h £ M h , 

(4.23) 

where 

/7i) = (acr™ 9 , or™ 9 ) + (“t(^t\ s - X h ), P$u™ 9 - y , h ) dT . 

(4.24) 


TeT„ 


Remark 4.5. Similar to the HDG methods, once A h is resolved, Uh and crh in (14.211) can he 
obtained in an element-by-element fashion. 

When applying the two-level algorithm, Algorithm [lj to WG methods based the model (14.231) . 
we set the operator Rh to be the symmetric Gauss-Seidel iteration or one sweep of Gauss-Seidel 
iteration. Similar to the HDG methods, one can easily show that the symmetric Gauss-Seidel 
iteration always satisfies Assumption IIIII 

When a is a piecewise constant matrix, from the HDG-like formulation (14.211) we can see that 
the WG framework (14.191) is essentially equivalent to the corresponding HDG framework. As a 
result, the convergence of the algorithm for he WG methods is as same as that for the corresponding 
HDG methods. 

For more general case of a , by using the same technique as in EDnacm it is easy to verify 
that 

a h {X h , A h ) ~ a™ 9 (\ h ,\ h ), VX h £ M h . (4.25) 

Then Assumptions [I] is obviously true for the WG methods. Following the same routines as in 
Section ET!?! one can derive the estimate (14.81) . Therefore, similar convergence results of Algorithm 
[Qfor HDG methods also hold for the WG methods. 

5 Numerical results 

In this section, we provide some numerical experiments in 2-dimensional case to support our 
theoretical analysis. We only consider Type 3 HDG method with olt = 1,VT £ Th ■ For more 
numerical results we refer to [la¬ 
in the first experiment, we set D = (—1,1) x (0,1)(J(0,1) x (—1,0] and define a(x,y) = 
diag (a(x,y),a(x,y)) with 

{ 1 , — 1 < x < 0 , 0 < y < 1 ; 

5, 0 < x < 1, 0 < y < 1; 

10, 0 < x < 1, — 1 < y < 0. 

Given an initial triangulation 7o of D, we produce a sequence of triangulations {7J : j = 1,2, • • • ,5} 

by a simple procedure: Tj+i is obtained by connecting the midpoints of each face of Tj for j = 
0,1, • • - ,4. 7o and 71 are presented in Figure |T] for clarity. For each 7) (j = 1,2, ,5), we 

set Th = Tj and construct TZh by using the standard V-cycle multigrid method based on the 
triangulations {% : i = 0,1, - - - ,j}, i.e. I — TZhAh denotes the error transfer operator of one 
V-cycle iteration. Here we set TZh and all smoothers encountered in the construction of TZh to be 
the symmetric Gauss-Seidel method with mo and mi iterations respectively. Using the standard 
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nodal basis for Mf ,, we let A/j be the stiffness matrix arising from the bilinear form (12.1211 . Suppose 
we are to solve Ahx = bh where bh is a zero vector, and we take xq = (1,1, • • • , 1)* to be the initial 
value, rather than the zero vector presented in Algorithm [T] We stop Algorithm Q] until the initial 
error, i.e. ^x^AhXo, is reduced by a factor of 10~ 8 . The corresponding numerical results (the 
number of iterations in Algorithm [[]) are presented in Table [1] 

The second experiment is a simple modification of the first one: we set IZh to be one sweep of 
Gauss-Seidel iteration. The corresponding numerical results are presented in Table [2] 




Figure 1: To (left) and 7i (right) 


Table 1: Numerical results for the first experiment 


k 

mo Ti 

r 2 

r 3 

r 4 

Ts 

k 

m 0 Ti 

r 2 

Ts 

Ti 

r 5 

k 

m 0 Ti 

r 2 

% 

T 

Ts 

0 

1 

19 

18 

19 

19 

19 

0 

1 

17 

18 

17 

17 

17 

0 

1 

17 

17 

17 

17 

17 


2 

13 

13 

14 

14 

15 


2 

12 

12 

12 

12 

12 


2 

12 

11 

11 

11 

11 


3 

10 

12 

13 

13 

14 


3 

10 

10 

10 

11 

11 


3 

10 

9 

10 

10 

10 

1 

1 

20 

21 

21 

20 

20 

1 

1 

20 

20 

20 

20 

19 

1 

1 

20 

20 

20 

19 

19 


2 

13 

14 

14 

15 

15 


2 

12 

13 

13 

13 

12 


2 

12 

13 

12 

12 

12 


3 

11 

12 

13 

13 

14 


3 

10 

10 

11 

11 

11 


3 

10 

10 

10 

10 

10 


mi = 1 mi = 2 mi = 3 


Table 2: Numerical results for the second exper¬ 
iment 


k 

mi 

Ti 

r 2 

r 3 

Ti 

Ts 

0 

i 

22 

24 

24 

23 

23 


2 

22 

23 

23 

23 

22 


3 

21 

23 

23 

22 

22 

1 

1 

34 

34 

34 

34 

34 


2 

34 

34 

34 

34 

34 


3 

34 

34 

34 

34 

34 


In the third experiment, we set = (—1, 1) x (—1, 1) and define a(x, y ) = diag(a(x, y), a(x , y)) 
with 

{ 1, — 1 < x < 0, — 1 < y < 0; 

7, 0 < x < 1, 1 < 2/ < 0; 

17, 0 < x < 1, 0 < y < 1; 

3, — 1 < x < 0, 0 < y < 1. 

We show the first two triangulations 7o and 71 in Figurc[5]and produce a sequence of triangulations 
{Tj : j = 0, ,25} in a successive way: Tj+i (j = 2,3, ,24) is obtained by refining the 
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smallest square containing the origin in Tj (in 71, the vertexes of the square to refine is in red 
color) as same as what has been done from 7o to 7i- T 25 is shown in Figure [3j The difference of 

the two-level algorithm between this experiment and the first one is that we simply take IZh = Ah 
here. The corresponding numerical results are presented in Table [5] 



Figure 2: To (left) and Tl (right) 



Figure 3: T 25 


Table 3: Numerical results for the third experi¬ 
ment 


k 

mo 

% 

Tio 

Tis 

720 

725 

0 

i 

15 

15 

15 

15 

15 


2 

12 

12 

12 

12 

12 


3 

12 

12 

12 

12 

12 

1 

1 

30 

30 

30 

30 

30 


2 

16 

16 

16 

16 

16 


3 

12 

12 

12 

12 

12 


For the first two examples, the regularity estimate holds with only a ^ 0.5, which violates 
the regularity requirement a € (0.5,1] in [19\ For the third example, not only (11.211 holds with 
a ^ 0.5, but also the triangulation is not quasi-uniform. However, for all the experiments, the 
numerical results are consistent with our theoretical results, which shows that our algorithm is 
convergent even when a is not greater than 0.5 in o and the triangulation is not quasi-uniform. 
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A Analysis of symmetric Gauss-Seidel iteration 

Let 7 Zh be the symmetric Gauss-Seidel iteration. As stated in Remark 13.31 we can show IZh 
satisfies Assumption Hill Suppose Assumption [I] is true. Then by the well-known properties of 
Gauss-Seidel iteration, we know that (13.131) holds with ut = 1. Thus it remains to verify (13.141) . 
Let {rji : i = 1, 2, • • • , N} be the standard nodal basis for M '/,. Define Pi : Mh spanj?^} by 

(-PjAft, Vi) Ah = {^h> Vi) Ah > * 1> 2, • • • , TV. (A.l) 

By Theorem 3 in [5], we have 


N 


{TZ h 1 ^h,^h)h = \\V'h\\ 2 Ah + inf 

E i= i !=1 


Pi ^ Mi 

j>i 


, 6 Mh- 


Ah 


Then, by using the same technique used in the proof of Lemma 3.7, we can obtain 

(Ph V-h,V-h)h 'y ] h T \\vh\\ h,dT ■ 

T&Th 

Denote 5ft := PhAh- By the definition dam of 7 Zh, it holds 

7 Zh = 272./J — PhAhPh = (25ft - 5 2 )Aft \ 


(A.2) 


(A.3) 


(A.4) 


which yields 

(7?.^ Aft, Aft)ft = (i?ft 1 S/ i (25ft — 5ft) 1 Aft,Aft)/ l . 

It is easy to verify that 5ft is symmetric with respect to the inner product (P~^ ■, •)/,. Then, from 
the inequality 

t(2t — t 2 ) -1 < 1 for all t S (0,1) (A.5) 

and the fact that all the eigenvalues of 5ft are in (0,1), it follows 

(Ph Aft, Aft)ft < (P^ 1 A/ l , Aft)ft, (A.6) 

which, together with (1A.3I) . leads to the desired result (13.141) . 
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